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We discuss the derivation of dissipative Bjorken hydrodynamics from a Schwarzschild black hole 
in asymptotically AdS spacetime of arbitrary dimension in the limit of large longitudinal proper 
time T. Using an appropriate slicing near the boundary, we calculate the Schwarzschild metric to 
next-to-next-to-leading order in the large r expansion as well as the dual stress-energy tensor on 
the boundary via holographic renormalization. At next-to-next-to-leading order, it is necessary to 
perturb the Schwarzschild metric in order to maintain boost invariance. The perturbation has a 
power law time dependence and leads to the same value of the ratio of viscosity to entropy density, 
l/(47r), as in the case of sinusoidal perturbations. Our results are in agreement with known time- 
' dependent asymptotic solutions of the Einstein equations in five dimensions. 
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I. INTRODUCTION 

There has been a substantial recent effort toward explaining the experimental results at the Relativistic Heavy Ion 

' ^ ' Collider (RHIC) in terms of a strongly coupled A/" = 4 Super Yang-Mills (SYM) gauge theory. While the plasma 

> formed at RHIC is described by a different gauge theory, namely Quantum Chromodynamics (QCD), one hopes that 

^ there are enough similarities between the two theories allowing us to aqcuire a good qualitative (if not quantitative) 

^ ' understanding of the behavior of the quark-gluon plasma. This will be achieved via the AdS/CFT correspondence 

, [3, Q which emerges in string theory and provides a method for studying the properties of strongly coupled A/" = 4 

■ SYM theory in terms of a dual gravitational description. 

' The AdS/CFT correspondence has been applied to the description of the quark-gluon plasma in several different 

^ • scenarios; see 0, |3| for nice reviews. Janik and Peschanski showed that a model for relativistic heavy ion collisions, 

00 \ suggested two decades previously by Bjorken ^5I|, can be seen to be a consequence of the AdS/CFT correspondence 

: They obtained an approximate solution to the Einstein equations in the five-dimensional bulk in the large longitudinal 

00 proper time (t) limit and showed that it gave rise to hydrodynamics on the four-dimensional boundary matching the 

' \ : Bjorken flow of an ideal fluid @. The work was subsequently continued to include subleading terms in the large t 

' expansion which could be understood as dissipative effects in the hydrodynamic behavior of the gauge theory plasma 

^\ 0,l,i[il. 

2? Kajantie, et al., showed that the two-dimensional Bjorken flow could be derived from a static black hole in 

the three-dimensional bulk spacetime by an appropriate time-dependent slicing near its boundary. We extended this 
result to arbitrary dimensions [ll| by showing that to leading order in r there exists a slicing near the boundary of an 
AdS Schwarzschild black hole which corresponds to Bjorken flow in the dual gauge theory plasma. The Schwarzschild 
metric to leading order in t agreed with the time-dependent asymptotic solution of Janik and Peschanski ^] in five 
' dimensions. In three dimensions it reduced to the result of Kajantie, et al., [T^ . 

The aim of this paper is to extend our earlier result [Jlj by including subleading corrections in the large r expansion. 
We show that next-to-leading-order corrections correspond to viscosity in the gauge theory plasma. At this level the 
coefficient of viscosity 77 is arbitrary, in agreement with results in five dimensions based on an asymptotic time- 
dependent solution to the Einstein equations At next-to-next-to-leading order we find that the Schwarzschild 
metric yields a flow which is not boost-invariant no matter how one chooses the slicing near the AdS boundary. Boost 
invariance is recovered after the Schwarzschild metric is perturbed by a power law r-dependent perturbation. We 
show that the perturbed metric is non-singular in the bulk, provided 

^ = - (1) 
s 47r 

where s is the entropy density, in agreement with asymptotic time-dependent solutions in five dimensions jSj]- This 
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special value of the ratio 77/s is also in agreement with the case of sinusoidal perturbations of an AdS Schwarzschild 
black hole [ij]. 

Our discussion starts with a review of dissipative Bjorken hydrodynamics in section[lll In section [TTTl we discuss the 
time dependent slicing we perform near the boundary in order to reproduce Bjorken hydrodynamics including next-to- 
next-to-leading-order contributions in the large t expansion. We also introduce the perturbation to the Schwarzschild 
metric which is necessary to maintain boost invariance at the order we are interested in. We show that demanding 
absence of singularities in the bulk metric leads to the standard value ([1]) of the viscosity to entropy density ratio. 
Finally section HVl contains our concluding remarks. 



II. DISSIPATIVE BJORKEN HYDRODYNAMICS 



Following Bjorken Q, in order to understand heavy ion collisions, we need to study boost invariant hydrodynamics. 
Let us consider a gauge theory fluid on a (d — l)-dimcnsional flat Minkowski space spanned by coordinates 
(/i = 0, 1, . . . , c? — 2). (We shall reserve the notation for the coordinates of a static gauge theory fluid; x^ will 
span the Minkowski space of the Bjorken fluid in order to avoid confusion.) With the colliding beams along the x^ 
direction, it is convenient to work with the coordinates r (longitudinal proper time) and y (rapidity) , defined by 



rcoshy , = rsinhy 

The {d — l)-dimensional Minkowski metric takes the form 



ds 



Bjorken 



= di^di^ = -dT^ + T^dy^ + [di^f 



(2) 



(3) 



where x-^ = (i^, . . . ^x'^ ^) represents the tranverse coordinates. 

For the stress-energy tensor we use the standard notions and enforce conservation and conformal invariance via: 



T^"" = (£+p)u^u'^+p.g^'^-CA''''VA'«^-r?(^A'^^VAu'' + A''^VAU^-^-^A'^''VAU^ 
V^T^'' = 

t; = 



(4) 



where /S^i, = g^^ -\- u^Ujy and e, p, r/ and C, represent the energy density, pressure, shear viscosity and bulk viscosity, 
respectively, of the fluid. Two constraints immediately follow from eq. 



{d-2)p , C = 



(5) 



In the rest frame of the conformal fluid, the velocity field is given by = (1,0). The stress-energy tensor simplifies 
to 









r) d-3 Vjr) 



V 

Choosing ly = t in the conservation equation of we obtain 



(6) 



d — I e ^d — 3 ri 

Assuming the viscosity to be a subleading effect, we deduce from ([7]) the leading behavior of the energy density 

^ ~ _(d-l)/(d-2) 



(7) 



(8) 



and hence of the pressure, temperature (from the Stefan-Boltzmann law in d — 1 dimensions, e ^ T'^ ^) and entropy 
density, respectively, 



P 



1 



d-2 d-2 Tid-i)/(d-2) ' 



To 



-l/(d-2) 



dp 
dT 



So 

T 



So 



1 £0 



d - 2 To 



(9) 
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The constants Eq and Tq represent the initial values of the energy density and temperature, respectively (at r = 1). 

At high temperatures we expect the viscosity to have the same dependence on the temperature as the entropy 
density (which is known to be true in the case of sinusoidal perturbations of the static A/" = 4 SYM plasma in five 
dimensions [l^) so that the ratio rj/s asymptotes to a constant. Therefore, we shall assume 

ry(r) « ^ (10) 

T 

where rjo is a constant. We may then solve eq. ([7]) and obtain a subleading correction to the energy density, 

^0 2770 



(11) 



r < 

yielding corresponding corrections to the temperature and entropy density, 



T ^ I 1 2r/o \ dp (I 2(d-2)r;o 1 



d-2 



Note that to leading order we obtain the ratio 

V _Vo _ {d - 2)-qoTQ 
s So (d-l)eo 



(13) 



This ratio is known to take the value l/(47r) (eq. ([T])) if the gauge theory fluid is dual to a perturbed AdS Schwarzschild 
black hole [l^. We shall show that the above Bjorken fluid also admits a gravity dual which is an appropriately 
perturbed AdS Schwarzschild black hole leading to the same value ([T]) of the viscosity to entropy density ratio. 

III. ADS SCHWARZSCHILD BLACK HOLE 

We are interested in obtaining the Bjorken flow described in section|lI]as a dual to a large AdS Schwarzschild black 
hole. The latter is a solution of the Einstein equations 

(^i?-Arf) (14) 

where A^j = _{!t:})^ti^l jg a negative cosmological constant. The solution best suited for our purposes is a large black 
hole with a flat horizon. Its metric can be written in the form 



1 - 2Aiz'i-i 

where x € W^~'^ and /it is an arbitrary integration constant. The Hawking temperature of the hole is given by 



(15) 



d-l 

= 4^ (1«) 

where is the horizon located at 

z+ - (2^^)-^ (17) 

and the boundary of the asymptotically AdS space is at z = 0. 

The VEV of the stress-energy tensor of the dual gauge theory can be constructed using holographic renormalization 
p^ . To this end, the metric needs to be brought in the form of a generally asymptotic AdS metric in Fefferman- 
Graham coordinates 

^^2 _ g^^dx^dx^ + dz^Q ^^^^ 

^FG 

Then near the boundary at zpc — 0, g^^i, may be expanded as 

5m. = 9^^} + 4g511^ + • • • + 4g + 0{zia) (19) 
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where = ry^^. For the black hole metric (|15p the transformation to Fefferman- Graham coordinates is achieved by 
writing the radial distance in the bulk as 



zfg 



^_rf_1^2/(rf-l) 



(20) 



[1 + 2^FG 



This brings the black hole metric in the form (|18|) with a diagonal g^i, where 

g^^^ ^^-l4~G? g^^^(l + !i,<l-}\''^'-'^ (z=l,...,d-2) (21) 

The dual gauge theory fluid is static and lives on the {d — l)-dimensional Minkowski space 

f^Sstatic = gj^dx^'dx" = rif.^dx'^dx'' = -df + dx"^ (22) 
Evidently, the first correction to the boundary metric is ^\ leading to the stress-energy tensor of the fluid 

where Gd is Newton's constant in the bulk. Explicitly, we obtain the energy density and pressure, respectively, 

.Hr")^.o, PHr-)^^, .o^^^ (24) 

obeying p ~ jZ2^i '<^^ expected for a conformal fluid (eq. (O). This is a perfect fluid and both coefficients of viscosity 
vanish (rj = ( — 0). With the temperature of the gauge theory fluid coinciding with the Hawking temperature (fTC)) . 
we deduce the equation of state 

1 ( A-kTh^ '^^^ 



and the energy and entropy densities, respectively, as functions of temperature, 

_ d-2 fAnTjA f^]^\ 
' 167rGd\d-lJ ' dTH~ AGd\d-lJ ^^^^ 

Since we are interested in Bjorken flow rather than a static fluid, we need to replace the transformation to Fefferman- 
Graham coordinates (PO)) with one that will yield a metric on the boundary {gj^J) of the form ^ rather than 
This is the task we turn to next. 



A. Leading order 



Near the boundary, the black hole metric ([15)) approaches 



dsli,. -> ^ (dsLtic + dz^) (27) 
leading to the boundary metric (1^^ . For Bjorken flow, we want instead, 

dsl.h. ^ (rfsljorken + di^) (28) 

leading to the boundary metric ([3]). We therefore need to find a transformation mapping the black hole coordinates 
(t, af, z) to a new coordinate system (r, ?/,a;',z) in a patch which includes the boundary relating the two asymptotic 
forms (f27|) and ([28|l . This will be done in the large r limit keeping the ratio 

^ = ZTJ^ (29) 
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fixed [y, [yj. More precisely, we define the large r limit as follows. Let 

r = To + r' (30) 

with constant tq ^ 1 and r' ~ 0{1), so that dr = dr' ^ 0{1). For the transverse dimensions, let ^ so 
that (dx^)'^ ^ C'(l)- The remaining term in the boundary metric ([3]) will be 0(1) provided we choose the rapidity 
y ~ C(1/t). Finally, for the bulk dimension let 



- l/(d-2) I -/ 



with zq, ~z' ~ 0(1) to ensure dz = dz' - 0(1) and v ~ 0(1). 
If we now apply the transformation [ll[ 



(31) 



d — 2 d--3 

t = -T''--^ 



d-3 



d — 3 
T-d-2, 



rl/(<i-2) 



(32) 



where a;-'" = (x^, ... 
z ^ 0), we arrive at 



„d-2 



), to the black hole metric pS]) (more presicely to a patch which includes the boundary 



ds 



b.h. 



(1 - 2^iv'^-^) dr^ + r^dy^ + {di^ f + - 



dz 



;2 



(33) 



with V defined in (|29p . The new coordinate system defines a new foliation near the boundary of z = const, hypersurfaces 
with time-dependent metrics leading to a flow on the boundary. At leading order in r, near the boundary the black 
hole metric behaves as in as desired. 

To see the form of the flow of the gauge theory fluid on the boundary, we need to bring (|33p in Fefferman- Graham 
form. This can be done to leading order in the large r expansion via the change of coordinates 



Z = ZpG 



'-FG 



d-l r{d-l)/{d-2) 



0{z 



2(d-l)x 
FG J 



(34) 



This brings the metric in the form ^5)1 with a diagonal g^i^ whose value at the boundary (g^i'J}) is given by ([3]). 
The first correction has non- vanishing components 



5^ = 



d-2 



1 



1 ^(d-l)/(d-2) 



(d-l) ^ 1 (d-l) ^ Ja^ ^ 

9yi ^2yyy ^_ ;l ^(d-i)/(d-2) 



(35) 



This metric agrees with the asymptotic time-dependent solution of the Einstein equations obtained by Janik and 
Peschanski f§\ in five dimensions. It leads to a stress-energy tensor for the conformal fluid on the boundary via 
holographic renormalization (eq. (j23p ) in agreement with expectations from Bjorken hydrodynamics (eqs. ([5]) and 
dH)), where 



{d~2)n 



(36) 



Notice also that the resulting energy and pressure agree with the case of a static boundary at initial time r = 1. 

The temperature of the gauge theory fluid can also be determined as it is related to the Hawking temperature (flB]) 
of the black hole. To see this, note that the static metric is conformally equivalent to the Bjorken metric 
Indeed, applying the restriction of the transformation (|32p on the boundary. 



we obtain at leading order in t. 



d — 2 d-3 
7" d-2 

d-3 



dsL 



T <i-2 y 



X 



^l/(d-2) 



2 

Bjorken 



-0(l/r)] 



(37) 



(38) 



The conformal factor, r ^/'^'^ scales the inverse of the Euclidean proper time period of the thermal Green function 
on the Bjorken boundary and hence the temperature of the gauge theory fluid, i.e.. 



T = 



-l/(d-2) 



(39) 
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where the proportionahty constant is chosen so that the temperature agrees with the temperature of the static fluid 
(Tff ) initiaUy (at t = 1). This is in agreement with the hydrodynamic result ^ with 



Tn = T, 



Finally, the entropy also agrees with ^ where 



So 



(2m) 



H 



(rf-2)/(d-l) 



4Gd 



The leading-order metric p3p can be seen to be regular in the bulk. Indeed, the Kretschmann scalar 

Tp2 _ n r/ABCD 

is found to be 



2(d-\) d + 2{d-2y{d-3)^i■'v 



2 2{d~l) 



.(d-3)/(<i-2) 



(40) 



(41) 



(42) 



(43) 



whose only singularity is obtained in the limit z — s- cx). 

Moreover, one can construct invariants which are linear combinations of the components of the Riemann tensor. 
For the geometry to be regular, these invariants must also be free of singularities [l7||. Introducing the vielbein 










( ' \ 

Z/T 




V ) 




z 


Vo/ 



(44) 



with a = 0, . 



, d — 1 spanning a local Minkowski space such that gAB^a^b ~ Vab- The Riemann tensor invariants are 



T^abcd = RABCDSaGf e^e'^ 



(45) 



Even though the individual components of the Riemann tensor have singularities, the invariants (|45[) are all regular. 
The non-vanishing components are found to be 



Ti-WW — "^2020 ~ ■ ■ ■ ~ 7^(d-2)0(d-2)0 = 1 + (d — 3)^1;''^"^ 
- ■ ■ ■ — 7?.(d-2)l(d-2)l = ■ ■ ■ ~ '^(d-3)(d-2)(d-3)(d-2) = -1 + ^ 
= • • • = 7?.(d-l)(d-2)(d-l)(d-2) = -1 - (d - 2i)liv'^^^ 

7^(<^-l)o(d-l)o = l-{d- 3)(d - 2)ixv''-^ 
together with those obtained using the symmetries of the Riemann tensor. 



Q(^ll^(d-:i)/(d-2)-^ 

e)(i/T-(rf-3)/(rf-2)^ 



(46) 



B. Next-to-leading order 

To extend the above results to next-to-leading order (©(l/r'^'*^'^^/'''^^^)), let us add a correction to the transfor- 
mation (I32p so that it reads 



d-3 2{d^2)T'^[l~2^lv<^-^)J T-{d-3)/{d-2) 

i_ I 1 Ci + bi{v) \ ±_.±f 1 Ci+ci{v) 

Ci 

A/(d-2) ^ 



(47) 



where Ci is an arbitrary constant, v is defined in (|29p and 61 (u), ci(u) and fi{v) are functions which vanish at the 
boundary (w = 0), so that they do not alter the boundary behavior of the metric obtained at leading order above. 
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The function fi{v) is determined by the requirement that the tz component of the metric vanish at the order we 
are interested in. We obtain the constraint 



w + (d - 2) (1 - 2nv'^-^y f[{v) = 
whose unique solution (with /i(0) = 0) may be written in terms of a hypcrgeometric function, 



-v^jd-S) 
2{d-2){d- 1 



-F I, 



d+1 



d- 1' 1 



2/111 



(d- 2)(d- 1)(1 - 2nv<i-^) 



(48) 



(49) 



With this choice of fi{v), under the transformation (|47l) the black hole metric (I15p turns into 



dsl.h. 



1 - 2fj,v' 
dp 



d-l 



d-1 



2{d~l)^j.Civ 

^(d-3)/(d-2) 



dr^ 



2biiv) 



^(d-3)/(d-2) 



T^dy^ 



2ci{v) 



r(d-3)/{d-2) 



1 _ 0,n,d-l I 2{d-l)tLCiV^ 



{dx^f 
(50) 



where the dots represent higher-order terms in the large r expansion. 

The Einstein equations at next-to-leading order yield two independent equations for the functions 6i(w), ci(w), 



b[{v) + {d - ^)c[{v) = 



-fiCiv 



d-2 , 



d-2 



[h{v) + (d - 3)ci{v)] + 



1 - 2fiv' 



d-i 



d-2 ' ' ' ' " {d-l){d-3) 
whose unique solution (with the boundary conditions 6i(0) = ci(0) = 0) is 



b[{v) 







bi{v) 



{d-3)Ci 



ln(l - 2fii 



Using (|52p. the metric (|50p can be written as 



ds 



b.h. 



1 + T-(<i-3)/(d-2) 

(di^)2 + 



ci(«) - ^ln(l-2/i«'*-i) 



dT^ + {l-2^,v^-'f-'^'''^' 



(51) 



(52) 



2 J„,2 



T 



+ (1 - 2/iw''-i) 



dP 

1 0,„..d-l I ■2{d-l)fj^Civ''-- 

± Z,/il7 -T (d-3)/(d-2) 



(53) 



which includes C'(l/r('^~'^)/^'^~^)) corrections to the leading-order expression p3|) . The next-to-leading-order expres- 
sion (j53p for the metric has no dependence on the rapidity y and transverse coordinates i"*", therefore it leads to a 
Bjorken flow for the gauge theory fluid on the boundary. 

We may now use holographic renormalization [l4| to calculate the VEV of the stress-energy tensor of the dual gauge 
theory. The transformation to Fefferman- Graham coordinates ([TSl) is 



Z = ZpG 



1-M 



Ci 



d-1 



r(d-3)/(d-2) J 



^d-1 
''FG 



^+0(z%-^)) 



(54) 



correcting the leading-order transformation (j34p . The form of the boundary metric is unaltered by design whereas 
the first non-vanishing correction away from the boundary reads 



2n{d~2) 



d-1 



1 



(d-l)Ci 



r(d-l)/{d-2) 



Jd-1) 



2/i 



1 



d-l T-{d~l)/(d-2) 



-,(d-l) _ 



1 



r2^VV 



2^ 

d-1 V r('*-i)/('^-2) 



{d-l){d-2)Ci 



(55) 



correcting the leading-order expression psp . 

Using eq. (|23p . we obtain a stress-energy tensor for the conformal fluid in agreement with Bjorken hydrodynamics 
(eqs. ([T0|) and pT|l ) with £o as before (eq. (p6|) ) and 



(rf- l)Ci£o 



(56) 
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matching the result of ref. Q for d = 5. 

The temperature of the gauge theory fluid can also be determined through the conformal factor relating the static 
metric (|22p to the Bjorken metric ([3]), as before. Applying the restriction of the transformation (j47|) on the boundary, 

t-J—^r^-^^C,\nr, ^=ry{-^^^-\, (57) 



we obtain at next-to-leading order in r, 

~ -l/(d-2) ^ ~r j L"*"Bjorkon 



^■^static — ( ^ ) ['^^Bjorkcn + ■ ■ •] (58) 



which yields the r-dependent temperature 



correcting the leading-order result (j39p and in agreement with the hydrodynamic result p2p with Tq = Th- The 
correct expression for the entropy density also follows and we obtain the ratio 

s Svr 

There is no constraint on this ratio at this order because the truncated metric ([55]) is regular in the bulk 'T]. This 
can be seen by a calculation of the Kretschmann scalar (|42)l . With the metric ((53l) . we obtain 



7^2 = 2(d- 1) 



d + 2(d-2)^(d-3)MV(-^) (l- J_3v|iS) ) 



(61) 



Eq. corrects the leading-order result showing that to this order the Kretschmann scalar is regular. 

However, the metric ([53)1 leads to singular Riemann invariants (|45|) . Indeed, we obtain explicitly, e.g., for d = 5, 

7^olOl = 1 + 2^,v^ + — ^ + . . . (62) 

exhibiting a simple pole at v = (2/i)^^/'*. This singularity should be absent, since our metric comes from a 
Schwarzschild black hole which has no singularities except as z — > oo. We obtained a pole because we have not 
included all contributions at order ©(l/r*^''"''^/^''"^^). There are additional contributions from next order {0{1/t)) 
terms are in the metric (I53p . Including them, the corrected metric reads 



"^b.h. — -2 



-[l-2^iv + _(,_3)/(,_2) ]dT +{l-2^iv ) T dy 



T 



r /^(d-3)/(ci-2) 



(1 - 2^1;'^-!) (d£^)2 + 1, ^ . 1 + 2yi^d£^dz + ^B^di^'dr + . . . '63) 

V ^ / ^ ' 1 _ 9,,„d-l I 2(rf-l)^Ciu'' ^ ' 

1 ^/j-u -r (d-3)/(£i-2) 



where the off-diagonal elements are 



, „ , [d-\)(d-Z)C^iirvv''-^ _ (d-l)CiM^^^"-^ 

' l-2//i;'i-i ' t(1-2^«^-i) ^ ^ 



These corrections do not lead to a Bjorken flow. However, the metric (|63|) satisfles the Einstein equations at this 
order. The Kretschmann scalar (|6ip is unaltered, and the Riemann invariants (|45p are corrected with the corrections 
cancelling all singularities. E.g., the invariant ([62]) for d = 5 is corrected to 



Teoioi = 1 + 2A.t'^ - -f . . . (65) 



which is a regular expression. 
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C. Next-to- next-to-leading order 



Extending the above results to next-to-next-to-leading order requires calculations which are considerably involved. 
We shall therefore restrict attention to the physically interesting case of five dimensions, setting d = 5, and employ 
Mathematica for the lengthy algebraic manipulations. The generalization to an arbitrary dimension is straightforward 
but adds little to the main results. 

Let us augment the transformation (|Tfl) for d = 5 with appropriate ©(l/r^/"^) terms as follows, 



t 



r2/3 



1 



r'^'y I 1 - 



t1/3 
V I 1 



Ci 



9(1-2^^4)^2 J ^1^1^^ 

Ci+b^jv) b2iv)+C2 
t2/3 r4/3 

Cl + Cljv) C2{v) +C2 
t2/3 r4/3 

02 (w) +C2' 



1^2 



r2/3 



/2(^) 
-4/3 



-2/3 



r4/3 



(66) 



where w is defined in The constant Ci is once again related to the viscosity coefficient (eq. ([55)) ). but to understand 
C2 one must employ second order hydrodynamics [15|. Whereby it is understood to be related to the relaxation time. 
The functions fi{v) and bi{v), ci{v) have already been determined at first perturbative order (eqs. p9)) and (|52p . 
respectively, with d — 5). The new functions f2{v) and 02 (w), ^(w), C2{v) ought to vanish at the boundary {v — 0) so 
as not to contribute to the boundary metric. 

As with fi{v), demanding that the tz component of the metric vanish at the order we are interested in yields the 
constraint on f2{v), 



3(1 - 2^iv'^Y f!2{v) - Ciw(3 -f IQ^iv^) = 



which has the unique solution (with /2(0) — 0) 

f2{v) = -Cih{v) + 



3(1 - 2^^4)2 



(67) 



(68) 



With this choice of f2{v), the application of the transformation (|66p to the black hole metric (fT5|) with d = 5 turns 
the latter into the form 



1 



- 1 - 2fiv-' 



1 - 



7-2/3 
2ci(«) , C2{v) 



ri/3 



dT^ 



2b,{v) , B2{v] 

7-2/3 



r4/3 



r'dy^ 



r2/3 



r4/3 



dz 



;2 



1 - 2^w4 



d2(v) 
7373" 



(69) 



where the dots represent higher-order terms. 
The off-diagonal elements are 



and we have defined 

A2{V) 

B2{v) 
C2{v) 



2,,2N 



3(1 - 2/ii;4)r4/3 



1 - 2fiv* 



r(l - 2fiv*) 



(70) 



_ ^^^4) - W(3C? + 2C2) - -(1 - 2^iv^)h{v) - 2(1 + 2^iv^)a2{v) 

bliv) - 2Cibi{v) - 202(1;) + 2b2{v) 
cl{v) - 2CiCi{v) - 2a2{v) + 2c2{v) 



d2{v) = 'i^iv*{3Cl + 2C2 + 2a2{v)) + 2v{l-2fiv'^)a'2iv) 



9(1 - 2^i;4) 



(71) 



Evidently, the metric depends on the rapidity as well as the transverse coordinates at next-to-next-to-leading order. 
This dependence cannot be eliminated by any choice of the functions which are yet to be determined. One may try 
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to modify the transformation ((66|) to eliminate the off-diagonal terms which depend on y and , but this only shifts 
the dependence on these coordinates to other components of the metric. If we insist on reproducing Bjorken flow on 
the boundary, we must perturb the Schwarzschild metric (fTSjl . Let the perturbed metric be 



ds 



perturbed 



— "*b.h. ~y 



T-4/3 



2Af^dx'^dz 



(72) 



where, apart from the off-diagonal elements, we are also modifying the zz component of the black hole metric by an 
amount proportional to an arbitrary function A{v). It turns out that, even though we have certain freedom in the 
choice oi A{v) (gauge freedom), this function cannot vanish. 

Using eqs. (j69p and (|72p. the zz component of the perturbed metric can be expanded as 



9~z 



\-2^lv^ r2/3(l-2/iw4)2 



v'D2{v) 

-4/3 



where 



D2{v) 



d2{v) 



A{v) 



(73) 



(74) 



v^{l-2^iv^Y (l-2/it;4)3 

Demanding that the perturbed metric satisfy the Einstein equations at the order we are interested in yields four 
constraints on the four functions ^2(f ), B2{v), C2{v), D2{v), 

3(1 - 2fivy{3 ~ 2^iv%B'^ + 2C^) - 3w(l - 2fiv'^f{B'^ + 2C^) 
-9v'{l - 2fxv^)'^D'2 + 18v(l + 6^iv'^){l - 2^iv^fD2 
+S^iv^[-l + mCl^iv' (ll-&^iv^ ~2{l-2^iv^)\n{l-2^iv'^))\ = 
9t;(l - 2^iv^fA'^ - 9(3 - lQpLV^){l - 2^iv'^fA'^ + 2%%^i'v\l - 2fj.v^)A2 
+9v'{3 - 2^j.v^){l - 2^iv'^fD'2 - 18i;(3 + %^iv^ - I2^''v^){l - 2^v^fD2 
18(3 + 2^u'*)(1-2^w4)3C2-8mw^ [7 + 2^w'* + 36CiW(9-f 44//v'' + 4//2t,8)] = o 

9w(l - 2^lv^fAl - 9(3 - 10/11''') (1 - 2iiv'^fA'2 + 2S^^?v'^ {I - 2fw'^)A2 
+9v{l - 2tiv^f{B'i + C'^) - 9(3 + 2nv^){l - 2nv^f{B'2 + C'^) 
-K9w2(3 - 2/Liw^)(l - 2y.v^ fD'2 - 18i;(3 + %[iv'^ - \2^'v^){\ - 2^v^ fD2 

~%^lv^{\^'\^[Lv^ ^-mcl[Lv'{m-2%^lv^ ^2%p?v^)\ = o 

72/ii;^(l - 2iiv^)A2 - 3(1 - 2[iv'^f{i - 2/iu^)(B^ + 2C2) - 36i;(l - 2[iv^fD2 

-f8/ii'^ [i'^ + 18Cf (l4/iz;4 + 4/iV - (1 - 2/it;4)(3- 2/iu4)ln(l - 2/ii;4))] =0 (75) 

coming from the rr, yy, xx and zz components of the Einstein equations, respectively. This system of equations does 
not completely determine the four functions. Keeping A2(v) arbitrary (gauge degree of freedom), the other three 
functions are determined to be 



-Hl8t;(l - 2iiv^fC'i 



-9(1 - 2^iv^f^ 



B'2{v) = 



A2{V) 

. + 2j.iv* 
2fiv^ 



9(1 - 4/i2^;8)2 

-72Ci^(l - 24:jiv 
+C3(1 - 2fiv*f - Ci{l - 2fiv*){3 + 4/i2t,8) 



[l - 2jw*){l + fiv* + 2fi'v'') 



20Ai\^) - 72C?(5 + 2fiv^ + 8ji''v^){l - 2jiv*) ln(l - 2/iw'*) 



2/1 



tanh ^ v' 



2/1 



A2{v) 
1 + 2/ti;4 

2jiv^ 

'9(l-4/l2-y8)2' 



- 2i;^(3 - ifiv* - 4/i^w«) - 2(1 - lOjiv" + 12n\f + Sji^v''') 



2/i 



tanh 



-1 2 

V 



-36Cf (11 - Gjiv* + 20ji\'^ + 2Aji'^v^^) - 36C^ (7 - 22jiv'^ + 20ii^v^ - Sfi'^v^') ln(l - 2jiv'^) 
+C3(1 - 2fiv^f + Ci{- '^ + 9fiv* - 10/t'i;' - 4/?„i2) I 



2/i 
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^2(f) = - 



'9(1 -2/Ltu4)(l + 2^^4)2 

2^z;2 

9(1- V„8)2(l_2^i;4) 



V?^2(l-6^ 



(1- V„8)2 

2 



tanh"^ v^y^i+ 108Ci ln(l - 2fiv'^) - ^ + C 



(76) 



The functions B2{v) and 6*2(1;) are found by integrating the first two equations, respectively. No arbitrary integration 
constants are introduced because -62(0) — 6*2(0) = 0. Notice also that apart from the arbitrary function A2{v), the 
above functions contain an arbitrary parameters C3 and C4. 

Constraints on the parameters are obtained by demanding regularity of the perturbed metric in the bulk. After 
some algebra, the Kretschmann scalar is found as an asymptotic expansion in t. 



7^2 



8{5 + 36^l'^v^) 

96^2^8 



2„,8 



230iCifi^v 
r2/3 



9(1 + 2/ii;4)T4/3 



IO8A2 + 



(2 - Sfiv'n 



3(1 - 2^iv^)C4 - 6(1 - 2^iv^) 



2fi 



tanh ^ 



(1 - 2^t-4)2 

2^+54Ciln(l-2^i;'') 



(77) 



correcting the lower order expression (|6ip for d = 5. At this order, we have a double pole at u = 1/ (2/i)^/4_ Demanding 
regularity of the Kretschmann scalar we obtain two constraints. One fixes the parameter Ci (which is related to the 
viscosity coefficient), 

= 6(^ ^^^^ 

and the other fixes the residue of the function A2{v) (which ought to have a simple pole at u = 1/ {2^)^^/^). Near the 
pole, we obtain 

Thus, A2{v) cannot vanish, however, other than the simple pole at w = l/(2/i)^/'*^ is arbitrary. Finally, there are 
no constraints on the parameters C3 and C4. 

The Riemann invariants (j45p are singular, as in lower order (eq. (j62p ). Their singularities are canceled by higher- 
order contributions to the metric (|72p which, however, spoil boost invariance (c/. with the corrected expression (j65p 
due to the corrected metric (1551) V 

Having obtained an explicit expression for the metric, we may now use holographic renormalization to compute the 
stress-energy tensor of the dual gauge theory. This is a tedious task. However, note that the temperature is easy to 
deduce from the restriction of the transformation (|66p to the boundary. 



t = -T^/^ 

2 



= T^'^y I 1 



2T2y2 _ (5^)2 



9r2 
Ci , C2 



— Ci In T — 



r2/3 



r4/3 



rl/3 



3C2 

2t2/3 



1 - 



r2/3 



-4/3 



(80) 



with a new correction dependent on C2 ■ The metric perturbation does not change the argument we employed at lower 
orders because the perturbation vanishes at the boundary. From the conformal factor relating the static and Bjorken 
metrics in the next-to-next-to-leading order we may obtain the temperature as 



T = T, 



H 



1 



Ci 



C2 



(81) 



The Stefan-Boltzmann law may be used to calculate the other thermodynamic quantities in the next-to-next-to-leading 
order. In particular, the viscosity to entropy density ratio is still given by (|60p with d = 5, but with Ci constrained 
by (|78p . It follows that this ratio is given by 
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which is the same value one obtains from sinusoidal perturbations of the AdS Schwarzschild metric [iSj . This is also in 
agreement with the conclusion reached by considering time-dependent asymptotic solutions of the Einstein equations 



IV. CONCLUSION 



We discussed the possibility of obtaining viscous Bjorken hydrodynamics on a (d— l)-dimensional Miskowski space 
from a large AdSd Schwarzschild black hole (of flat horizon). The latter is normally considered dual to a static 
gauge theory fluid on the boudary whose temperature coincides with the Hawking temperature. By appropriately 
modifying the boundary conditions, we obtained viscous Bjorken hydrodynamics on the boundary in the limit of large 
longitudinal proper time (r — > oo) at next-to-leading order. Our results are in agreement with those obtained by 
considering time-dependent asymptotic solutions of the Einstein equations in five dimensions Moreover, since 

our bulk space consisted of a Schwarzschild black hole, we were able to determine the temperature of the conformal 
fluid on the boundary in terms of the Hawking temperature of the hole. At next-to-next-to-leading order, we saw 
that no choice of boundary conditions could lead to a boost-invariant flow. In order to obtain a dual Bjorken flow 
at that order, we had to perturb the black hole metric. This led to a constraint on the viscosity coefficient and the 
viscosity to entropy density ratio was fixed to the value l/(47r) as in the case of sinusoidal perturbations flS*!. This 
was in agreement with a next-to-next-to- leading order calculation of a time-dependent asymptotic solution ^] . 

We should point out that a constraint on the viscosity is not necessary if one does not perturb the black hole metric. 
In this case, one obtains deviations from Bjorken flow which are a subleading effect at late times. It might be worth 
exploring the connection of such deviations (coming from a dual Schwarzschild black hole) to existing (RHIC) and 
forthcoming (LHC) experimental data. 

It may also be worthwhile, albeit tedious, to go beyond the perturbative order considered here. It has already 
been observed that the supergravity Fefferman-Graham metrics dual to boost-invariant hydrodynamics suffer from 
singularities of the curvature invariants near the reputed black hole horizon [T6| . This occurs at the third perturbative 
order it the large r expansion and cannot be canceled by an appropriate choice of transport coefficients as has been 
done at second order. However, by working with Eddington-Finkelstein instead of Fefferman-Graham coordinates, an 
expansion in the new time coordinate was shown to lead to non-singular solutions at all orders with the correct choice 
of transport coefficients (itL . It would be interesting to investigate the connection of the AdS Schwarzschild black 
hole metric with these Eddington-Finkelstein-type solutions of the Einstein equations. 
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